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SUMMARY 


The present paper is an extension of a previous investigation 
(described in NACA TN 2383) concerned with the solution of the nonlinear 
differential equation for transonic flow past a wavy wall. In the 
present work several new notions are introduced which permit the solu- 
tion of the recursion formulas arising from the method of integration 
in series . In addition, a novel numerical test of convergence, applied 
to the power series (in the transonic similarity parameter) representing 
the local Mach number distribution at the boundary, indicates that smooth 
symmetrical potential flow past the wavy wall is no longer possible once 
the critical value of the stream Mach number has been exceeded. 


INTRODUCTION 


In NACA TN 2383 (ref. l) the writer considered the problem of two- 
dimensional transonic flow past an infinitely long sinusoidal solid 
boundary. The problem was treated in the physical plane and the purpose 
was to investigate whether or not the flow past the wavy wall remains a 
smooth symmetrical type of potential flow when the undisturbed -stream 
Mach number exceeds its critical value. By a smooth type of potential 
flow is meant one for which the velocity potential, say, and its first 
derivatives are single-valued and continuous; that is, there are no 
discontinuities of the nature of shock waves. 

Initially, the Prandtl-Busemann small-perturbation method was 
applied and the velocity potential developed inclusive of the third 

power in the "thickness" coefficient e = f 


where 

a 


X 


amplitude of wavy wall 
wave length of wavy wall 
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The velocity potential was then referred to the critical speed of sound 

kfl - m2) 3/2 

c cr , the coefficient € was replaced by ILJ. — s_ where 

7 + 1 

(7+l)€ 

k = is the transonic similarity parameter, and terms 

involving powers of 1 - Me 2 higher than the first were neglected 

(main assumption for transonic flow) . This simplified or transonic 
form of the Prandtl-Busemann solution was shown to be identical with 
the one obtained by solving the simplified nonlinear differential equa- 
tion for transonic flow, with the boundary condition taken not at the 
wave-shaped wall but at the flat plate corresponding to vanishing 
amplitude. The calculation was carried through the sixth power in the 
transonic similarity parameter k and corresponds to the insuperable 
task of obtaining the Prandtl-Busemann solution to the sixth power in 
the thickness coefficient g- Thus each iteration step of the Prandtl- 
Busemann method contributes to the transonic form of the solution, which 
may therefore be considered a result of thin-profile theory with disturb 
ances not necessarily small compared with 1 - Moo 2 . The main conclusion 
reached in reference 1 was that the transonic similarity parameter k 
must be less than 4/3 — a value still somewhat greater than the critical 
value 0.83770 calculated there. 

The purpose of the present work is to express the solution of the 
problem of transonic flow past the wavy wall in a form more suitable for 
general considerations and to prove that the assumed smooth symmetrical 
type of potential flow cannot exist at stream Mach numbers beyond the 
critical value. 


SYMBOIS 

x, y nondimens ional rectangular Cartesian coordinates 

a amplitude of wavy wall 

X, wave length of wavy wall 

X./2n reference element of length 

e "thickness" coefficient of wavy wall, 

7 ratio of specific heats at constant pressure and constant 

volume 

U stream velocity 
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M local Mach number 

Moo undisturbed-stream Mach number 


k 

c cr 

0 

f n(y) 


transonic similarity parameter, — — 

critical speed of sound 
velocity potential of flow 
transonic disturbance potential 
functions of y only, related to f(x,y) 


numerical coefficients 


generating functions of k, E a" P L 


2r 


Ajjj p functions of dummy variable r, A^_£ r n } lower label 

n=l,m 

starts from 1 when m is negative and from m when m is 
positive 

Primes denote differentiation with respect to independent variable. 


ANALYSIS 

General Formulas 


If the undisturbed stream is in the direction of the positive x-axis, 
then the velocity potential 0 referred to the critical velocity c cr 
can be written as (see ref. 1) 


7Tl(l - M„ 2 )f(x,y) 


0 


X + 



k 


NACA TN ZjkQ 


where the second term on the right-hand side is a disturbance velocity- 
potential and implies that terms involving powers of 1 - higher 

than the first have been neglected. The differential equation for 
transonic flow satisfied by the function f(x,y) is obtained from the 
general differential equation for compressible flow and takes the 
following simplified nonlinear form: 


d 2 f _ bf d£f 
dy 2 c)x 2 


The boundary conditions to be fulfilled by f(x,y) are as follows: 


■\ 



(at y = 00) 


r ( 2 ) 


bf 

by 


-k sin x 


(at y = 0 , -00 < x < <») 


Here x and y are nondimens ional rectangular Cartesian coordinates 
simply related to the physical plane coordinates X and Y by means of 
the transformation 


x = X 


y = (1 - M oa 2 ) 1 ^ 2 Y 


and the equation of the wavy wall is Y = € cos X or 

7 - (l - ^ e cos x. Clearly, f (x,y) involves only the variables 
x and y and the transonic similarity parameter k. 


The most general form for the function f(x,y) to ensure symmetrical 
potential flow past the wavy wall is the following (see fig. l) : 
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f<x,y) 


-X 



f n sin nx 


(3) 


where the f n are functions of y only. When this form for f(x,y) 
is substituted into equation (l) and the coefficients of the separate 
harmonic terms sin nx are equated to zero,, the following system of 
second-order nonlinear ordinary differential equations for f n results: 


n-1 


- n 2 f = - 
n 


1 


IZ m ( n " m ) f m f h-m - | n H “U + m)f, 

m=l ^ m=l 


m^n+m 


(n - 1, 2 , . . . oo) (4) 

Before proceeding to the solution of these equations, several formulas 
of general interest and subsequent use are given. They have been derived 
in reference 1. The local Mach number distribution is given by 

; 

1 - M 2 = -(l - Moo 2 )g ( 5 ) 


The equation from which the critical value of the transonic simi- 
larity parameter is calculated follows from equation (5) by taking 
x = 0, y = 0, and M = 1; that is, 



= 0 


or with the aid of equation (3)^ 


00 

H nf n(°) = 1 


n=l 


( 6 ) 


The pressure coefficient 


n _ P ~ 

^p,M„ I 72 
' 00 ^Poo U 
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is given by 


c P,M„=-rrr( 1 - M » 2 )( 1 + |) (7) 


Integration in Series for the Functions f 

In reference 1, equations (4) were solved by an iteration pro- 
cedure and approximate expressions for f^ to fg were obtained. An 

examination of these expressions showed that the general form of f n is 

n-1 

oo 2p+n -2 00 

-(2p+n)y Y - qT np n+2r 

f n = Z_ e Z_ y A q r k 

p=0 q=0 r=p 

(n = 1, 2, . . . 00 ) (8) 

where, if p = 0, the upper limit of q is n-1 and, if p ^ 0, the 
upper limit of q is 2p + n - 2. The four-labeled coefficients A^ P 

are real numbers calculated from recursion formulas obtained from the 
system of differential equations (4) and the boundary condition at the 
surface of the wavy wall. The boundary condition at y = 00 is auto- 
matically satisfied by the form of whereas the boundary condition 

at the wall takes the form 


.\ -k 

(n = 1) 

*N 

Cr 

n 

0 

M 

0 

(n ^ 1) 



(9) 


Inserting the expression for f n given by equation (8) into equations (9) 
yields immediately the following results: 



= 1 


Y~ (2p + n)A^ P 
L 0 r 

-n-n 


r 

E 

p=o, 1 


A? p 

1 r 


r do) 


(n = 1, 2, ... 00 ) 

J 
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where, if n = 1, the lower limit of p on the right-hand side is unity 
and, if n ^ 1, the lower limit of p is zero. Also, if n = 1, the 
upper limit r of p goes from 1 to «> and, if n ^ 1, r goes 
from 0 to 00 . 


By elementary 
can be replaced by 


manipulation of series, the second of equations 
the following more useful forms: 


( 10 ) 


nA n 0 = A n 0 
0 0 10 


(n = 2, 3, • 


nA 


n 0 
0 r 


A 1 r - L (2p + n)A 0 r 
p=0, 1 p=l 


(r = 1, 2, 

\ n = 1, 2, 



(11) 


where in the first term on the right-hand side the lower label p starts 
from 1 when n = 1 and from 0 when n / 1. 


In reference 1, recursion formulas were derived for the coef- 
ficients A^ p . In the present paper a much more significant approach 

is introduced. Note that equation (8) can be rewritten in the following 
form: 


= O 

p=0 


ke 


_y\2p+n 


n-1 

2p+n-2 

r 

q=0 


A V ^ 


where 


( 12 ) 


A n P 

q 


r A n p k 2r 

r=0 q P+r 


In a manner similar to that described in reference 1 for the coef- 
ficients A^ p , recursion formulas can be obtained for the power 
n p 

series A . Indeed, the two types of recursion formulas are intimately 

connected and, for a given value of p, the one can be obtained from the 
other by mere inspection. A single recursion formula can be written for 
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the general quantity A n q ^ hut the resulting expression is cumbersome 

and serves no practical purpose. It is much more desirable to obtain 
separate recursion formulas for each value of p. As examples of pro- 
cedure, the recursion formulas for p = 0, 1, and 2 are considered in 
the sections ' which follow. 


Recursion Formula for A n ^ 

With p = 0, the recursion formula is (compare with eq. (57) of 
ref. 1) 


_ , , , „n 0 , , w „>„ii-2 n 0 

2n( q + l)A q+1 = (q + l)(q +2)6 q A q+2 + 


where 



q n ' 2 - q l 

^ ^ (m + l)(n 

Q 1 =0 m=q-q 1 

(n = 2j 3j • • • 00 } 


- m - l)A n “ m_1 0 A m+1 0 
- q-Qi 

q = 0 , 1 , . . . n - 2) (13) 


n-2 0 (q = n - 2) 

5 q = 1 (q / n - 2) . 

This recursion formula can be solved, the solution starting with 
q = n - 2 and descending towards q = 0. Thus, for q = n - 2, equa- 
tion ( 13 ) becomes 


n-2 

8(n - 1 )aJJ_° = r (m+ 1) (n - m - l)A 
m=0 


m+1 0 ^n-m-1 0 
m n-m-2 


(n = 2, 3, • .••.“) (14) 


Now, multiply both sides of this equation by r 11 where r is a dummy 
variable and sum from n = 2 to n = oo. Then 


n-2 


8 ^ (n - l)A*_° r n = ^ Y1 (m + ' 1)(n ' m " ^^m ° A *" 111 " 1 °’ 


n-m-2 


n=l 


n=2 m=0 


(15) 
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Let 


A 1 0 = XI A n-l rI1 


n=l 


then 


00 

A 1 O’ = L <-°i 
n=l 


and 


00 n-1 

( A 1 0*) 2 = / 1,11-2 V m ( n - m)A m 0 A n_m 0 

i—- m-1 n-m-1 


n=2 m=l 


By observation, it can be seen that the right-hand side of equation (15) 
is equal to 0 *)^* Equation (15) can thus be replaced by the 

following first-order nonlinear differential equation: 


( rA l o ') 2 - 8rA l o' + 8a i 0 = 0 


(16) 


The solution of this equation is 


A 


1 0 


c o e 


-p A I 

A i 0 


where c Q is the arbitrary constant of integration. From the defini- 
tion of A-]_ q it follows that, with r = 0, 


c = A 1 0 
c o A 0 


ancT hence 


J 


A 1 0 

A 1 0 
A 0 


= e 


*L rA l 0 Al °* 
? Oio 
A 0 


(17) 
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Note that by definition A n ® involves the coefficients Aj? ^ which 

are ultimately calculated by means of the auxiliary relations, equa- 
tions (11), engendered by the boundary condition at the surface of the 
wavy wall. In solutions of the recursion formulas, therefore, the 

coefficients A 11 ^ appear as undetermined quantities. 

A 1 o' 10 

The expansion of in powers of rA . is a nontrivial 

a t u 

A 0 

problem which fortunately can be solved with the aid of Lagrange's 
investigations on the reversion of power series. In reference 2, 
Lagrange's problem is illustrated by the following example: 

To expand e ax in powers of y = xe . The result is expressed as 


^ax _ 


= 1 + ay 


a(a - 2b) o + a ( a ~ 3b)‘ 
2* y 3! 


+ 


a (a - rib) 11 " 1 n 


n: 


(18) 


In particular, with a = 1, b = -1, and % = e x , the solution of the 
equation log | = y| or £ = (generally referred to as Eisenstein's 

problem) is given by the series 


i 


00 

L 

n=l 


n 


n-2 


n-1 


(n - 1)! 


Then, if £ is replaced by 
of equation (17) is 


"1 0 
il 0 


and 


y by 


1 rA 1 0 

If 0 > 


the solution 


00 
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f- nA" 0 r"' 1 = f- l£± f A l 0\” r n 

L- n_1 U ( n - lJjV^X 0 ' 


Hence, by equating coefficients of equal powers of r on both sides of 
this equation. 


n"~ g / A 1 O')" 

..lx"- 1 ' 0 ) 


(n - 2, 3, . . . oo) 


which is the solution of the recursion formula, equation (l4). The 
n 0 

coefficient A n _ 1 can be considered as the generating function for the 
set of coefficients A n _^ r . Thus, equation (19) can be written as 


a n 0 ? 2r _ n 11 ” 2 V s - i 0 2r\ , 

Z_ A n-1 r k --TH=l 2_ A 0r k (20 

r =0 nlk J 

Then, by equating coefficients of equal powers of k on both sides of 

this equation and putting aJ ° = 1 (see boundary condition, eq. (10)), 
the following equations are obtained: 

A " 0 n n-2 1 

"- 1 0 _ 

n!4 


n 0 = . n 11 " 2 1 0 

Vl 1 " (n - DU "' 1 0 1 


A 11 0 

^n-l 2 


(n - 1) 


- A 1 0 + _ n n ~ 2 / i 0 \£ 

!4 n_1 0 2 2!(n - 2)1^1° V 
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2 0 1 

Note that from the first of these formulas, A = — and that from 

10 8 

the boundary conditions, equations (11), a q 0 = l£' 

Consider now q = n - 3* Then equation (13) becomes 


n 0 


2n(n - 2 )aJJ_ 2 = (n - l)(n - 2)A^_J + | n ^ m(n - m)AjJ_° (22) 

m=2 

\ 

Multiply both sides of this equation by r 11 and sum from n = 3 to 
n = 00 . Then 


n-1 


2 H n(n - 2)A^_° r n = T (n - l)(n - 2 )a“_° r n 
n=2 n=l 


Let 


5 II nr ' 1 fl m(n - m)A I-2 CI-1 (23 > 

n=3 m=2 


*1 0 


A n ° r 11 
n-1 r 


and 


Ao.-f A n0 r n 
A 2 0 ~ L A n-2 r 


n=2 


then 


00 n 


n=l m=l 




A 1 o' A 2 o' 
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It can be shown easily that the last term on the right-hand side of 
equation (23) is equal to | r^A-j. Q 'k 2 0 *^ * . Hence, equation (23) 
can be rewritten as the following differential equation: 


2r 2 A 2 0 " - 2rA2 0 « = 0 " - 2xk ± Q « + 2A X Q + | Q 'A 2 0 ')' 


Now, from equation (l6), it follows that 


(24) 


( rA i o ’) 2 " 8 ( rA l o’ " A 1 0) 


hc)'=a 


4A n 0 * 


rA l o' 


rA 


1 0 


>'( rA l o') ” 

Then equation (24) can be written as 


= 4rA " 

4rA l 0 


(25) 




6rA i o' + ‘‘hj 


(26) 


Let rA 2 0 ' = v 


K o' )* 


4 - rAi 0 

Then equation (2 6 ) becomes 


-, where v is the new dependent variable. 


•( rA i o') 2 ^ 1 * 1 °'^ ’ 6rAl °’ + tAl 



Ik 
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or with the aid of equations (25), 


Hence 


" ■ IK o')' 


v = — r A 


. . 1 c 

8 1 0 + 8 1 


or 


(rA-, n *)3 (rA, *\ 2 

8rAp n ' = ^ T + c-i f — 1 P * , 
20 4 - rA-j_ q' 1 4 - rA^ q' 


The arbritary constant c^ is determined by differentiating this last 
equation twice and evaluating for r = 0. Thus, 


a 2 o° 

c n = 64 0 


K°) : 


Finally, 


\2 0 


'2 o' “15 r2 [[ A l 0') 2 ]’ + 8 —^2 rA l o" 

r o j 


(27) 


Now, from equation (l8) with a = 2, b = -1, and .1 = e x , it follows 
that 


I = 2 


n' 


n-3 


n-2 


(n - 2)! 


n=2 
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Or with 


A ' 
1 0 

.1 0 


1.10 

y = n rA 0 


n=2 


and 




Thus, equation (27) becomes 


I 1 *- 2 

n=2 


n-1 


00 



n 


n-3 


(n - 3) !4 


n-2 



+ 


8 



L - i)a s-? 


n=2 


r 


n-1 


n 0 

Hence, replacing a“_ 1 by its value from equation ( 19) and equating 

the coefficients of equal powers of r on both sides of this equation 
yields 


*2-i = 8 



+ 


1 

2 


n 


n-k 


(n - 3) !4 : 


n-1 



n 


(n = 3, b, . . .00) (28) 
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the solution of the recursion formula, equation (22). Again, equa- 
tion (28) can be written as 



Then, equating the coefficients of equal powers of k on both sides of 

1 0 
0 0 


this equation and putting A^ ® = 1 and ^ yields: 


. n 0 
l n-2 0 


n 0 
l n-2 1 


n 0 
n-2 2 


n 


n-3 


rr 


n-4 


2(n - 2) !4 n_1 2(n - 3) !4 n ~ 1 


8n n " 3 


A^? + 


n 


n-3 


,1 0 


, . ,.n-l 0 1' n-1 “0 1 

(n - 2) !4 (n - 3) !4 


8n 3 T A§ 2 + — t A n p + 8A n ? A n 

(n -2)!4 n - 1 02 (n- SJ^L 0 2 0 1 ° 


> (30) 


(n - 3) 
n "3 LI 0\2 


1/ A 1 of + ^ (a 1 o\ 

2 ^ 01 ^J 2(n - 4) !4 n-1 


Note that from the first of these formulas A p 0 = 'Tff and that from 
the boundary conditions, equations (ll), a q 0 = ^>* 
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The procedure followed in order to obtain equation ( 28 ) for A r ^ 
is a general one and with very little difficulty other members of the 
family A n ^ can be obtained. 

For example. 



From this equation it follows that 


,n-4 


A n 0 _ 2 n^ + 1 

n ' 3 U 3 ( n _ 3 ) * 4 n_1 2 ( n -4)'4 n_1 


, n 0 |,o n n-lf a 3 0 _ D n 11- ^ „2 0 . S. / -> n \ 

A n-3 1 = 48 ^TI A 0 1 + 8 Aq 1 + | (32) 


(n - 3)14* ' (n - 4) ! 4* 

I 1 + 1 L_J 

[3 (n - 4) ! 2 (n - 5 ) LJ 


n 11- ^ 10 

, n -1 0 1 


N 


From the first of these formulas and the boundary conditions, equa- 
tions ( 11 ), aJ ° = 3 ^ and A^ ° = 

At this point, it is noted that the coefficients of the form A^ q 

are calculated from the first formula of each set, equations ( 21 ), ( 30 ), 
(32), and so forth. A number of this type of coefficient have been y 
evaluated (see ref. l) by means of the recursion formulas for the 
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coefficients A a ^ 


themselves. They are listed as follows: 


2 0 
1 0 


a3 0 
1 0 


4 o 
1 0 


1 

A 6 0 - 

7 x 13 

8 

A 1 0 " 

72 x 256 

1 

a 7 0 _ 

13 

2? 

A 10’ 

18 x 256 

. 7 

A 8 0 - 

13 X 19 

384 

A 1 0 " 

576 x 256 



7 

768 


A careful examination of these numerical values leads to the general 
rule. 


nO l 3n - g} 

10 n-1 

n!4 


(n — 2, 3> ... °o) 


and from the boundary conditions, equations (ll). 


V n 0 
0 0 


{3n - 5} 


nn 


Ik' 


n-1 


( n — 2^ 3> * • • 00 ) 


where by definition 

{3n - 5} = 1 X 4 x 7 X 10 x 13 X . . . x (3n - 5) 

In the expression for the local Mach number distribution evaluated 
at the crest of the wavy wall (x =0, y = 0), there occurs the following 
power series: 

00 

F - V nA§ § k n 
n=l 
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or 


F 



This power series can be expressed in the closed form 


F = 2 - 2 


i 1 -x*> 


\2/3 


The graph of F against k is a semicubical parabola with the cusp 
k 

point at k = — and F = 2. With the necessary condition that one and 

only one value of k correspond to a given value of F, the transonic 
similarity parameter k cannot be greater than 4/3- Moreover, the 
lower limit of k is zero when the amplitude of the wavy wall is zero 
but Moo is different from unity. 


Recursion Formula for A 


n 1 

q 

,n 1 


With p = 1, the recursion formula for A n is (compare with 
eq. (69) of ref. l) 

Mn + 1)AV - 2(n + 2)(q + 1)6^ A^ + (q + l)(q + 2)S n -J' n Aj + | = 


_q_ n-l-q 1 
q 2=0 m=q-q 1 


q q+1 

n-m 0ml 1 


- | nS q YL IZ " (n " m)An "”l A ”-^l ' 2 n(n + 1 )a1 0° A " +1 ° 


where 


(n = 1, 2, . . . ooj q = 0, 1, . . . n) (33) 


R n-l,n 0 

6 9 1 


( q = n - 1 or n) 
( q ^ n - 1 or n) 


K n ° 

5 q " 1 


(q = n) 
(q / n) 
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The solution of this recursion formula proceeds as in the case p = 0/ 
the solution starting with q = n and descending towards q = 0. Thus, 
for q = n, equation (33) becomes 





A n+1 0 


or 


(n = 1, 2, . • . 00 ) 


A 


n-1 1 
n-1 



DA 1 0 A n ° 
1;A 0 n-1 


(n = 2, 3> • • • 00 ) 


.n 0 


Hence, inserting the expression for A n _ 1 given by equation (19) gives 

,1 0\ n+1 


A- V - - I -AlfT (n - 2, 3, . . . -) (3M 

n_1 8 (n - 2)\k n - 1 \ ° ' 


or 


1 = . i ° n ~ 3 l Y 2 ° k 2r \ 

n - lr+1 s (n . 0r / 


n+1 


From this relation, by equating coefficients of equal powers of k on 
both sides of the equation, the following equations are obtained 

.n-1 1 _ _ 1 nI1 ~ 3 


n_1 1 8 (n - 2) !U n_1 


, n-1 1 


n 


n-3 


,10 1 n 1 


n-3 


,1 0 


n- 1 2 8 (n - 2) ’A 11-1 01 8 (n - ° ± 

> r, n_ 3 


1 n-1 1 _ 
l n-l 3 " 


! ( A o i) 2 ] - h 


D- 0 / A i 0 ^ 

2l _ 1 

n 11 - 8 T 

-|A 0 2 + (A q ly ) 

J 8 

(n - 3) !4 n_1 


°\ 2 


n - 4) ’.i)- 11 * 1 ° 

1/ 

( n = 2, 3, . 


0 


> (35) 



NACA TN 27 b 8 


21 


Consider now q = n - 1 > equation ( 33 ) becomes 


n-1 


4 (n+l)A“ } = 2 n(n + 2 )A n 1 - inL (m + 1 ) (n - m - l)A^ v A 
n ~ ± n d m=0 


m+1 0 .n-m-1 1 
m n-m-1 


1 / n x . 1 0 . n +1 0 

- n(n + l)A 0 n _ 1 


{n = 1 , 2 , . . . ») 


Multiply both sides of this equation by r n and sum from n = 1 to 
n = °°. Then 


4 V” (n + 1)AJJ j r 11 = 2 n(n + 2)A n n 1 ^ " 

n=I ’ n=l 


oo n-1 

1 's n V” / , w xx.m+l 0 .n-m-1 1 

- Z_ nr 11 l_ (m + l)(n - m - l)A A 


2 n=l m=0 


m n-m-1 


00 

1 A 1 ® n(n + l)A n+ \° r n 

2 0 n-1 


( 36 ) 


Let 


'0 


1 * 1 1 r “ 


n=l 


A -Y~ A n 1 T n 
A 11 L A n-1 r 
n=l 


V - > A n 0 r 11 
10'L n-1 r 

n=l 


A 


2 0 


= Ti A “ 2 pn 

n=2 


Then it can easily be shown that the second term on the right-hand side 
of equation (36) is - | r^r 2 A Q 1 , A 1 ' and that equation (36) can 

be replaced by 


/ 



<?> 
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H rA l l)' = l" , + 6rA o i' - 5 r { T \ i’ A i o')' 


- rA 1 0 A " 
2 0 2 0 


Now, 


*0 1 = - 5 a 1 0 ° £ <" - - 11 - 1 


n=l 


or, with the aid of equations (25)> 


(37) 


rA = _ Jl. a 1 °frA ' 'V 2 

rA o l 64 0 \ 1 0 / 


Hence, from the use of equations (25) and equation (27) for A 2 0 ' it 

follows, after routine calcinations, that the integral of equation (37) 
is 


rA n 


i = w a1 o°(- i( rA i o') 2 - *( rA i o') 3 - ^r( rA i o') k + 
( C 1 + 4 ) jK rA i o') 2 ■ 2r \ o]} 


where 


C;L = 64 


* 2 0 


(cy 


Now from equation (l8), 

( 


m 


rA ' ] = m 

1 0 


V n n - m - 1 f A l 01“ 
4j ( n - 0 1 


„n 


(38) 
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Then 



Equating coefficients of equal powers of r on both sides of this 
equation gives 



(n = 2, 3, • • • °°) (40) 

From this equation the expressions for all the coefficients of the 
type . with n = 2, 3> • • • 00 and r = 0, 1, . . . oo can 

be obtained. Thus, for example. 
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n-1 1 1 n n ~^r 20 . 70 . 37 ^ 4 ~ | 

n-2 1 " 32 4 n_ l li. n - 2)1 (n - 3) ! (n - 4) ! (n - 5)]J 


A n-1 1 = _ A. SL-l JfT^n + Da 1 ° + 32A 2 4 

n-2 2 32 4 n-l | ±)A 0 1 ii2A 0 lj(n - 2) ! 


|7jn + 1 )aJ ° + 32A§ ^ (n -^ 3) - , + [(37n + 5)Aj ° 


r (4D 


256a 2 0 


0 1 ,/«10p.2o\ 1 

lj(n - 4)! + V^O 1 + 8A 0 lj(n"'5). 




For q = n - 2, equation (33) becomes 


4(n + 1)A“_1 = 2(n - l)(n + 2 )a£_ 1 - n(n - l)A n ^ 
n-2 _ 


d m=0 


/ ~w „\ .m+2 0 . n-m-2 1 

(m + 2)(n - m - 2)A m A n _ m _ 2 + 


(■ - D(» - - - UA“i 0 A”--y] - 


1 , lh 10 A n+1 0 

2 n(n + 1)A 0 A n-2 


(42) 


The solution of this recursion formula follows along the same lines as 
that for equation ( 36 ) and leads to the following result: 
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n ? 

Recursion Formula for A 

With p = 2. the recursion formula for A n „^ is 
' q 

8(n + 2)A n 2 - 2(n + 4)(q + l)6 n+2 A n f + (q + l) ( q + 2)8 n+1,n+2 A?J 
q q q+1 q q + ^ 

n +l-q 1 


1 a n+2 E n-1, n, n+1 -0,1,2 r 9- T — , ^ , 

-2 n8 q 5 ’q 6 m Z_ Z_ ‘ 2 >< n " "■ + 

q 1= 0 m=q-q x 


2)A n '“ +2 0 A m-2 2 _ i nS n+l,n+2 8 n y- m („ . „) A ' A - 

q l q ' q l 11 q q lqy0m=7q 1 

\ n(n + l)(s n+2 A 1 . 0 A n+1 1 + 6 n+1 ’ n+2 A 1 . 1 A n+1 0 + 

2 \q°q q Oq 


.n-m 1 A m 1 


6 °,n + 2 A 1 1 A n + 1 0\ . n(n + 2) U+2 2 0 A n + 2 0 & 0 A 2 0 A n + 2 o\ 

q 1 q-l / \ q 0 q q 1 q -1 j 


(n — 1, 2, • . • 00 ^ q — 0, 1, • • • n + 2) ( ^5) 


where delta is defined in the usual manner. The solution of this 
recursion formula proceeds as in the previous cases, starting with 
q = n + 2 and descending towards q = 0. Thus, for q = n + 2, equa- 
tion (45) becomes: 


or 


»n 2 
n+2 


nA £ 


0 A n+2 0 
L n+1 


(n = 1, 2, 


A n-2 21/ 2)a 2 0 n 0 

R n 8^ ’ 1 n-1 


(n = 3, 4, 
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Then, from equation ( 19 ), 

A n-2 2 _ 1 , 

6 b K 


n 


2 ) _^L/ a 1 0\ 

' n!4“-A A 0 1 


n+2 


or 


*n-2 2 v 2r _ (n - 2)n : 

A n r + 2 k Tn+2 

n !4 


(n = 3, 4 , 


n ‘ 2/ ^ A 1 0 k ^" +2 
A 0 r K 



By equating coefficients of equal powers of k on both sides of 
equation, the following equations are obtained: 


A n " 2 

2 _ 

(n - 2)n n " 2 

n 

2 

i ), n+2 

n :4 

A n-2 

2 

(n 2 - 4 )n n " 2 

n 

3 

i i, n+2 

n !4 

A n-2 

n 

2 

4 " 

(n 2 - 4 )n n-2 

, ( ,n+2 

n !4 


Consider now q = n + 1} equation ( 45 ) becomes 
8(n + 2 )a£ + 2 = 2(n + 2)(n + 4 )a£ + 2 - 


n-2 


1 ^ / , lW x m+1 0 n-m-1 2 

5 n 2 _ (m + l)(n - m - 1 )A m A n _ m+1 


m=0 


1 / -i \ / * 1 0 * n+1 1.11. n+1 o\ 

— n( n + 1 ) f A o A n+1 + A ± A n J - 

N r 

( , 0 \ / a 2 0 a n+2 0 , a 2 0 A n+2 o\ 

n(n + 2 )(A Q A n+1 + A x A n j 

(n = 1, 2, . . . 


4 ( 46 ) 


this 


( 47 ) 


oo) ( 48 ) 
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With the introduction of 


, A n 2 r n 

l—, A n+0 r 


A 


-2 2 tn n+2 
n=l 


= A n 2 r 11 
-12 L n+1 r 
n=l 

_ V" A n 0 r 11 

1 0 2 _ n-1 

n=l 


00 

n 0 n 

A 2 0 = z_ A n-2 r 


n=2 


A = A n 1 r 11 
A 0 1 l— R n r 
n=l 


equation (48) is replaced by the following ordinary linear differential 
equation: 

8 ^r^A_ 2 _ 2 ) = 2r 2 (rA _ 2 2 ') ' + 12r^rA_ 2 2 'j + l6rA_ 2 2 - 

1 

1 —2/2. in t\ _ i 0 r 2n " _ i ^ r^A " - 

2 V 1 0 " 2 2 ) 2 0 r A 0 1 2 1 10 


A 2 0 _n " . A 2 0 » | _ U r A " + 

■ 0 rA l 0 A 0 A 10 A 1 rA 2 0 


v2 0 


20 a , ,1 0 ,2 0 


A c u A » _ A A 
1 2 0 0 0 


(49) 


By repeated use of equations (25), the solution of this differential 
equation is found to be as follows: 


8A -i a - ( Al o°) 2 [i( rA i o') - js{ TA i o ') 2 - ^a( rA i o') 

aSraK o ') 4 + T& A i o'*5 # 2o* - *) [|( rA i o') : 

r2j 4 o]| - Y r ( Al o°) + a 2 o° fio' g( rA i of]- 


rA 1 0 A 2 0 
0 0 


(50) 
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Finally, 

A n ' 2 2 = 
A n-1 . 


with the aid of 



5fe( c ! + V- 


equation (38) and the definition of A_j_ 2> 

1 2 ~ ) n n -3 _ n n ~5 

( n - 2)! (n - 3)|j4n-l ~ 64 ( n _ i^.^n-l + 




From this equation, there follows in the usual manner formulas for the 
coefficients of the type A^ 2 Thus, for r = 0, 


A n-2 2 
n-1 2 



1 2 ~ j n n ~8 n 11- ^ 

( n " 2 ) '• ( n " 3)^Jj^n+3 _ 4) i h n+ 2 

1 ~ 1 n 11 " 2 

(n - 2)jJJn?2 = 3 ’ ' 


. 00 ) 


DISCUSSION OF CONVERGENCE OF SMOOTH TYPE OF POTENTIAL 
FLOW PAST WAVY WALL 


In the preceding sections, a number of examples of recursion 
formulas and their solutions have been given in considerable detail. 

The purpose of this exposition is threefold: First, to show the inherent 

elegance of the method of integration in series although the equations 
concerned are nonlinear in character; second, to present a type of 
analysis which may be useful in other problems involving nonlinear 
differential equations; and third, to indicate that an analytical proof 
of convergence may ultimately be obtained by careful examination of the 

recursion formulas for the quantities A n g p and their solutions. One 


example is the obtaining of the general 
.n 0 

A^ and the subsequent conclusion that 


expression for the 



coefficients 


In actual practice it has been found more convenient to evaluate 
the coefficients A^ p from their separate recursion formulas rather 
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than to derive the general formulas. The appendix contains the exact 
numerical values of the coefficients necessary for the development of 
the functions f n to the eighth power in the transonic similarity 

parameter k. These coefficients are utilized to demonstrate numerically 
the test for convergence of the smooth symmetrical type of potential flow 
(eq. (3)) assumed in this paper. Thus, when the form of f n given by 
equation (8) is inserted into equation (5) for the local Mach number 
distribution and evaluated at the surface of the sinusoidal wall (y = 0), 
the following result is obtained: 


M 2 - 1 
1 - Moo 2 


= -1 + y 2 " k n S (n - 2m)cos(n - 2m) x \ ^n-2m p (52) 

x <- — 0 m 


n=l m=0 


where denotes the integral part of n/2. At the crest of the wavy 

wall (x = 0), the point of maximum fluid velocity, with the numerical 
values for the coefficients inserted, equation (52) becomes: 


+ ! = k + 1 k 2 + k 3 + J2L + k 5 


1 - V 


576 X 256 


359381 6 . 7326757 v 7 A 81688733 , 8 


k + . . . 


270 x 256 


6480 x 256 


864oo x 256 


The critical value of k (that is, when M = 1) calculated from this 
equation is 

k cr = 0.83244 


(Note that in ref. 1 the value . kc r = O.83770 corresponds to the first 
six terms of eq. (53)-) Consider now the infinite series 


00 

T. A" ? 
n=l 00 


where. 


{3” - 5} 


(n = 2, 3> • • • «.) 
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The Cauchy ratio test R ln = ■ ■ ^ ^ yields in the limit n— the 

A 0 0 

result that the radius of convergence is equal to or less than 4/3. 

If the corresponding ratios R 2 n are formed for the right-hand side of 

Rl 

equation (53) and the quotient N n = = — is calculated, the resulting 

K 2n 

sequence of numbers is as follows: 


n 

^ R ln 

R 2n 

N n 

1 

16.0 

8.0 

2.0 

2 

4.5 

1.92 

2.34375 

3 

3.04762 

1.78041 

1.7H75 

4 ' 

2.5 

1.33366 

1.87472 

5 

2.21546 

1.35000 

1.64104 

6 

2.o4l67 

1.17721 

1.73432 

7 

1.92481 

1.19588 

1.60950 


The noteworthy feature of this table is that although R^ n (and pre- 
sumably R2n) converging quite slowly toward R^ = — (and R 2 ), the 

quotient N n exhibits a strong tendency to approach an asymptotic 

value N for a relatively low value of n. Figure 2 shows this tendency 
in a graphic manner. The apparent asymptote represented by the straight 
line is the ratio of 4 / 3 , the limit of Rj n as n— ♦ and of 0 . 83244 ,' 

the critical value of k. Certainly, the rapid approach of the lower 
dotted curve toward the apparent asymptote and the decreasing oscilla- 
tions represented by the upper dotted curve indicate that the critical 
value of k is the radius of convergence R 2 of the power series on 
the right-hand side of equation (53). Thus, the critical value of the 
stream Mach number marks the limit of convergence of the smooth sym- 
metrical type of potential flow assumed. The ability to approx ima te 
closely the limiting value N is a matter of luck} namely, the choice 
of the known comparison series. Once, however, the proper comparison 
series has been selected and the approach to an asymptote indicated, there 
can be no question of the meaningfulness of the approximate value of N 
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obtained. It may be that one would like to extend figure 2 to n = 9* 
This extension would entail the forbidding calculation of an additional 

185 coefficients A^ The result presumably would be to. decrease 


slightly the approximate critical value of k and thereby raise slightly 
the straight-line asumptote of figure 2. This extension of figure 2 
would show still more convincingly the approach to an asymptote and the 
dying-out of the oscillations. Perhaps more important still, figure 2 
definitely shows that conclusions based on less than six or eight terms 
are mere speculations in this field. 


CONCLUDING REMARKS 


If the numerical test of convergence presented is acceptable, the 
conclusion to be drawn is that smooth symmetrical potential flow past 
the wavy wall exists only for the purely subsonic range. Moreover no 
such flow can possibly represent the transonic or mixed type for which 
a local region of supersonic flow near the solid boundary is imbedded 
in the otherwise subsonic stream. It follows as a corollary that the 
transonic or mixed type of flow past the wavy wall is necessarily an 
asymmetric one. This asymmetry in the flow pattern entails a resistance 
usually defined as wave drag. As shown by experimental observations, the 
shock wave associated with wave drag closes the downstream portion of 
the local supersonic zone. 

As a final remark — although the analysis and conclusions of the 
present work refer directly to the wavy wall, the suggested result that 
the critical stream Mach number marks the limit of smooth potential flow 
very likely applies to other boundaries. This conclusion is based on 
the idea that the gradual transition from a purely sinusoidal wall to a 
boundary composed of a single bump, say, introduces no essential changes 
in the analysis. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 6, 1952. 
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APPENDIX 


NUMERICAL VALUES OF THE COEFFICIENTS A n p 

q r 


Coefficients for f^: 


A 


1 o _ 

oo" 


i 



A 1 2 = 65 

02 576 x 256 


A 1 3 = 338? 

0 3 1152 x 256 2 


A 1 0 - JL 

0 1 256 


0 2 


2765 

9 x 256 s 


A 1 2 = 26435 

03 81 x 256 3 


A 1 3 = 3907 

13 432 X 256 2 


A 10, 1861 

0 2 " 3 X 2 5 6 2 


Ai 1 - - 

v 3 " 


184345 
1728 x 256 2 


A 1 2 = 
1 2 


7 

72" x 256 


A 1 3 = 23 

23 12 x 256 2 


A 10 s 4896755 

03 81 x 256 3 


11 l 

L 1 1 = " m 


1 1 = 33 

12 " 6^ x 256 


A 1 2 = _ 53995 

1 3 1728 x 256 2 


A = - Z 
2 2 32 X 256 


A 1 3 = . 583 

33 54 x 256 2 


A f 3 a . m 

k 3 12 x 256 2 


A 1 1 - . 1.39, A 1 2 - - - A 1 3 - 1 

1 3 " 1, x 256 2 • 2 3 " - 96 x 256 2 ^ 3 " " 96 x 256 


A 


1 2 _ 
3 2" 


1 

8 x 256 


A 1 2 = - ^ 

3 3 8 x 256 2 
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Coefficients for f 2 : 


A 


2 0 = 1 
oo 16 


A 


2 1 _ 125 

01 " 36 x 256 


.22. 122U5 A 2 3 _ , 17999 

*° 2 144 x 256 2 ^ 3 90 x 2563 


»2 0 _ jj!9 
01 72 X 256 


.21 = . 45895 

02 108 x 256 2 


a 2 2 = 35023583 
0 3 1620 x 2563 


A 2 3 

A 1 3 


72577 

5760 x 256 2 


a 2 0 = 234215 

02 432 x 256 2 


a 2 0 = 42791533 
03 648 x 256 3 


2 1 = 4210409 

0 3 = ‘ 81 x 256 3 


A 21 a 5 

A 1 1 256 


a 22 = 6245 

12 72 X 256 2 


a 2 2 = 104 q63 
13 90 x 256 3 


a 2 3 e 45409 
2 3 2160 x 256 s 


a 2 3 _ 3515 

33 864 x 256 2 


0 

ft l 0 


1 

8 


21. 815 

1 2 ' ' 576 X 256 


a 22 = . 12 a? 3 s _ Ji .0322 

22 3 x 256 2 ^ 3 864 x 256 2 


A 


2 0 


1 1 


11 

1024 


2 i = _ 456185 A 2 2 = - 1 6 7 

1 3 3456 x 256 2 2 3 4 x 256 2 


A 2 3 = _ 73 

53 8 x 256 2 


.2 0. 4085 .2 1 1 

1 2 ' 24 x 256 2 2 1 ' " 128 


a 2 2 _ _ 47 a 2 3 = . 2 3 

3 2 192 x 256 63 12 X 256 2 


A 2 0 = 21287 A 2 1 - - 11 

13 324 X 256 2 2 2 ~ 32 x 256 



24199 

864 x 256 2 


2 1 _ _ 4811 2 2 = _ l 

23 192 x 256 2 k 2 12 x 256 


2 

L 4 3 ° ' 


11 

2 x 256 s 
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Coefficients for f^: 


. 3 o 

1 

a 3 l 

1765 

A 3 2 - 

475823 

A o 0 = 

72 

A 0 l " 

3 x 256 2 

A 0 2 " 

15 X 256 3 

ft 3 o 

23603 

A 3 1 - 

1351393 

A 3 2 - 

74353 

A o 1 = 

27 x 256 2 

A 0 2 

2880 x 256 2 

A 1 2 

320 X 256 2 

.3 o 

9843883 

A 3 ! = 
1 1 

155 

A 3 2 = 

3083 

A 0 2 = 

81 x 256 3 

32 x 256 

A 2 2 

24 X 256 2 

A 3 0 = 
1 0 

1 

A 3 1 

30365 

A 3 2 - 

4l 

2S 

A 1 2 - 

48 x 256 2 

A 3 2 - 

2 x 256 2 

.3 0 

259 

A 3 1 = 
A 2 1 

117 

A 3 2 = 

1785 

A 1 1 = 

72 x 256 

32 x 256 

4 2 

4o x 256 2 

.3 0 

149317 

A 3 1 - 

29545 

A 3 2 = 

25 

A 1 2 = 

432 x 256 2 

A 2 2 

96 x 256 2 

5 2 

2 x 256 2 

A 3 0 = 
2 0 . 

l 

32 

A 3 1 = 
3 1 

l 

" 256 



A 3 0 - 

33 

A 3 1 = 

. _55_ 



A 2 1 ' 

32 x 256 

3 2 

256 s 




.3 0 _ 139 

22 2 X 256 2 
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Coefficients for f^: 


A 4 0 = 

7 

A 4 1 = 

-35251 

• a 4 2 _ 

741115 

0 0 

6 x 256 

A 0 1 

90 x 256 2 

A 0 2 " 

5184 x 256 2 

II 

O rH 

^<P 

3905^7 

a 4i = 
0 2. 

80872738 

A 4 2 - 
1 2 " 

38381086 

720 x 256 2 

675 x 2568 

405 , x 2563 

11 

O OJ 

1704729613 

A 4 1 = 

2465 

A 4 2 - 

423833 

16200 x 256 3 

1 1 

576 x 256 

2 2“ 

1152 x 256 2 

a 4o = 

7 

A 4 1 = 

30596538 

A 4 2 - 
3 2 

20845 

1 0 

384 

1 2 

135 x 256 s 

144 x 256 2 

A 4 0 = 

- 32947 

> 

ro 4 =~ 

H M 
II 

" 455 

A 4 2 

173 

1 1 

36 x 256 s 

96 x 256 

A 4 2 = 

12 X 256 2 

A 4 0 = 

20527213 

A 4 1 - 

127825 

A 4 2 

• ‘ 30 

1 2 

162 x 256 3 

A 2 2 

192 x 256 2 

A,_ _ = 
5 2 

256 2 

A 4 0 = 

3 

A 4 1 = 

119 

A 4 2 - 

36 

2 0 

128 

3 1 

48 x 256 

6 2 

5 x 256 2 

A 4 0 = 

617 

A 4 1 - 

51343 



2 1 

288 x 256 

A 3 2 - 

216 x 256 2 



A 4 0 = 

369587 

A 4 1 - 

25 



2 2 

1728 x 256 s 


48 x 256 



A 4 0 = 

l 

A 4 1 - 

275 



3 0 

96 

a 4 2 ' 

8 x 256 2 



A 4 0 = 

11 





3 l 

24 X 256 





A 4 0 - 

4811 




: ' . 

3 2 

144 x 256 s 







NACA TN 27 48 


37 


Coefficients for f^: 


A 5 0 = 1 
0 0 ' 3840 


a 5 o = 90354659' 

0 1 259200 x 256 2 


A 5 

A 1 


0 = 7 
o 7S8 


A 5 1 _ _ 5470394 

01 8l x 256 3 

a 5 l = . 88735 

11 96 x 256 2 

A 5 1 _ _ 25 x 15587 

21 288 x 256 2 


A 5 

A 1 


0 

1 


886243 

1080 X 256 2 


A 5 1 _ _ 25 x 997 
3 1 24 x 256 2 


5 0 _ _25_ 
2 0 1536 


A 5 1 = _ -105 
4 l 6 h x 256 


a 5 0 = 24375 

21 32 x 256 s 


A 5 1 = _ 9 

51 32 x 256 


A 5 ° _ _5_ 

3 0 384 

a 5 0 = 8950 

31 27 x 256 s 


A 


5 

4 


0 = 25 

0 24 x 256 


A 5 0 = 1375 ' 

^ 1 24 x 256 2 
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Coefficients 

for f£: 



6 0 

91 

A 6 1 - 

2483629 

A o 0 

432 x 256 

A o 1 

54 x 256 3 

a 6 o 

285522407 

A 6 1 = 

4371541 

0 1 

4860 x 256 3 

1 1 

5760 x 256 2 

a 6 o 

91 

a 6 1 _ 

265901. 

A i 0 

72 x 256 

A 2 1 " 

192 x 256 s 

. 6 0 

144862403 

A 6 1 - 

403535 

A i i 

810 x 256 3 

A 3 1 “ 

288 x 256 s 

A 6 0 

23 

A 6 1 - 

• 79841 

ii 

2 0 

8 x 256 


96 x 256 2 

A 6 0 
A 2 x 

2553091 
2880 x 256 2 

A 6 1 - 
5 1 

2205 

8. x 256 s 

A 6 0 

_ 13 

A 6 1 = 

2401 

3 o 

1024 

6 l 

60 x 256 2 

A 6 0 

56387 



3 1 

96 x 256 2 



A? 0 
4 0 

15 

8 x 256 



A 6 0 
\ 1 

815 

4 x 256 2 



A 6 ° 
5 0 

9 

20 x 256 



A 6 0 

297 



A 5 1 

10 x 256 2 
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Coefficients for fy: 


A T 0 _ 13 

00 126 x 256 


a 7 0 = 13 

10 18 x 256 


a 7 0 = 287 

2 0 HM X 256 


A 7 0 = 833 

30 288 X 256 


A 7 0 = 3^3 

^ 0 ikk x 256 


A 7 0 = 3^3 

50 320 x 256 


a 7 0 _ 2^01 

60 o 

45 X 256 
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Coefficients for fg: 


,8 0 = 13 x 19 
00 18 x 256 2 


a 8 0 = 13 X 19 
1 0 576 x 256 


A 8 0 = 89 

2 0 64 x 256 


A 8 0 = 59 

3 0 24 x 256 


8 0 _ *7 

4 0 18 x 256 


A 8 0 =-J£- 
5 0 2880 


A 8 0 = 7 

6 0 2880 


A 8 0 - _1_ 

7 0 2520 
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